MAN/HEP/2009/5 



OS 
O 
O 

(N 
C 

>—> 

O 

(N 



a,; 

CD 



> 

o 

o 

cn 

o 

On 
O 



X 



Central Exclusive Production at the LHC 

J.R. Forshaw a * 

a School of Physics & Astronomy, University of Manchester, Oxford Road, Manchester M13 9PL, UK. 

After a brief resume of the theory underpinning the central exclusive process (CEP) pp —* p + H + p, attention 
is focussed upon Higgs bosons produced in the Standard Model, the MSSM and the NMSSM. In all cases, CEP 
adds significantly to the physics potential of the LHC and in some scenarios it may be crucial. 



1. Introduction 

The idea to install detectors far from the inter- 
action point at CMS and/or ATLAS with the ca- 
pacity to detect protons scattered through small 
angles has gained a great deal of attention in re- 
cent years and the report presented in [T] consti- 
tutes a significant milestone on the road to CEP 
physics at the LHC. fn this talk, I should like 
to focus attention in particular on Higgs boson 
production, as illustrated in Fig. [TJ For a much 
more extensive survey of the physics that can be 
studied after installing forward detectors see pQ. 
Detection of both protons by detectors located 
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Higgs 
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Figure 1. Central Exclusive Production of a Higgs 
boson. 



420m from the IP at ATLAS/CMS has the virtue 
that central systems with masses up to 200 GeV 
can be observed with an event-by-event precision 
of around 2 to 3 GeV, and adding also 220m de- 
tectors extends the reach to much higher masses. 
The clean environment of CEP generally makes 
for reduced backgrounds (even in the presence of 



significant amounts of pile-up) and that is often 
aided by the fact that the centrally produced sys- 
tem is predominantly in a J z = 0, C-even, P-even 
state. Of course having such a spin-parity filter 
also provides an excellent handle on the nature of 
any new physics. For very little extra cost, for- 
ward detectors promise to significantly enhance 
the physics potential of the LHC. 

That said, CEP is not without its challenges. 
The theory is difficult, triggering can be tricky, 
signal rates for new physics are often low, new 
detectors need building and installing, and pile- 
up needs to be brought under control. With the 
arrival of data on CEP from CDF at the Tevatron, 
confidence is building in the theoretical modelling 
and extensive studies have demonstrated that all 
of the other challenges can be met, e.g. see pQ. 



*Talk presented at the workshop "New Trends in HERA 
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Figure 2. The relevant lowest order Feynman di- 
agram for qq — > q + H + q. 



Without further ado, let us very quickly re- 



1 



2 



J.R Forshaw 



view the Durham model for CEP, more details 
can be found in |2|3j . The calculation starts 
from the easier to compute parton level process 
qq — > q + H + q shown in Fig. [5] The Higgs is 
produced via a top quark loop and a minimum 
of two gluons need to be exchanged in order that 
no colour be transferred between the incoming 
and outgoing quarks. The real part of the am- 
plitude is small and the imaginary part can be 
determined by considering only the cut diagram 
in Fig. [2j The calculation of the amplitude is 
straightforward: 



d(PS) 2 6(( qi ~Q) 2 )6((q 2 +Q) 2 ) 



2gq? 2 ff(Z2Q 2gtf 2gq$ 



Q 2 



k\ 



^2 



yab c c a b 
' [lv im'jn'mk'nl 



(1) 



We write Q — aqi+Pq 2 +Qr- The delta functions 
hx the cut quark lines to be on-shell, which means 
that a « -(3 » Q r /s < 1 and Q 2 « Q\ = -Q^. 
As always, we neglect terms that are energy sup- 
pressed such as the product a/3. In the Standard 
Model, the Higgs production vertex is 
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where V = m 2 H a s j '(4irv)F(m 2 H /m 2 ) and F w 2/3 
provided the Higgs is not too heavy. The Durham 
group also include a NLO K-factor correction to 
this vertex. 

We can compute the contraction q\V^q 2 ei- 
ther directly or by utilising gauge invariance, 
which requires that k%V$ = k%V$ = 0. Writ- 
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system should have a vanishing z-component of 
angular momentum in the limit that the protons 
scatter through zero angle (i.e. q%^Q\)- Since 
we are interested in very small angle scattering 
this selection rule is effective. One immediate 
consequence is that the Higgs decay to 6-quarks 
may now be viable. This is because, for massless 
quarks, the lowest order qq background vanishes 
identically (it does not vanish at NLO). The lead- 
ing order bb background is therefore suppressed 
by a factor ~ m 2 /m 2 H . 

Returning to the task in hand, we can write (y 
is the rapidity of the Higgs) 
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We are mainly interested in the forward scatter- 
ing limit whence 
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As it stands, the integral over Qt diverges. Let 
us not worry about that for now and instead turn 
our attention to how to convert this parton level 
cross-section into the hadron level cross-section 
we need. 



x, Q T -fc r 




dG(x,Q T ) 
d ln(Q T ! ) 



since 2k\ ■ k 2 rs x\X 2 s ~ m-H- Note that it is as 
if the gluons which fuse to produce the Higgs are 
transversely polarized, e; ~ k^. Moreover, in the 
limiting case that the outgoing quarks carry no 
transverse momentum Qt = —k±T — k 2 T and so 
ei = — e 2 . This is an important result; it general- 
izes to the statement that the centrally produced 



1 We can do this because Xi 
Sudakov components are ~ ( 



mnlsjs whilst the other 



Figure 3. The recipe for replacing the quark line 
(left) by a proton line (right). 



What we really want is the hadronic matrix ele- 
ment which represents the coupling of two gluons 
into a proton, and this is really an off-diagonal 
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parton distribution function. At present we don't 
have much knowledge of these distributions, how- 
ever we do know the diagonal gluon distribution 
function. Fig. [3] illustrates the Durham prescrip- 
tion for coupling the two gluons into a proton 
rather than a quark. The factor K would equal 
unity if x' = x and kr = 0, which is the diagonal 
limit. That we should, in the amplitude, replace 
a factor of a s Cjr/7r by dG(x, Qt) / dh\Q T can be 
easily derived starting from the DGLAP equa- 
tion for evolution off an initial quark distribution 
given by q(x) = 6(1 — x). The Durham approach 
makes use of a result derived in [4] which states 
that in the case x' <€. x and k T <C Q\ the off- 
diagonality can be approximated by a multiplica- 
tive factor, K . Assuming a Gaussian form factor 
suppression for the fcr-dependence they estimate 
that 
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and this result is obtained assuming a simple 
power-law behaviour of the gluon density, i.e. 
G(x,Q) ~ x~ x . For the production of a 120 
GeV Higgs boson at the LHC, K ~ 1.2 x e-^ll 2 . 
In the cross-section, the off-diagonality therefore 
provides an enhancement of (1.2) 4 « 2. Clearly 
the current lack of knowledge of the off-diagonal 
gluon is one source of uncertainty in the calcu- 
lation. The slope parameter b is fixed by assum- 
ing the same fc^-dependence as for diffractive J jt\> 
production^, i.e. b ps 4 GeV -2 . 

Thus, after integrating over the transverse mo- 
menta of the scattered protons we have 
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where f(x,Q) = dG(x,Q)/d\nQ 2 and we have 
neglected the exchanged transverse momentum in 
the integrand. 

Now it is time to worry about the fact that our 
integral diverges in the infra-red. Fortunately I 
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have missed some crucial physics. The lowest or- 
der diagram is not enough, virtual graphs possess 
logarithms in the ratio Qt/^h which are very 
important as Qt — > 0; these logarithms need to 
be summed to all orders. This is Sudakov physics: 
thinking in terms of real emissions we must be 
sure to forbid real emissions into the final state. 
Let's worry about real gluon emission off the two 
gluons which fuse to make the Higgs. The emis- 
sion probability for a single gluon is (assuming for 
the moment a fixed coupling a s ) 
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The integration limits are kinematic except for 
the lower limit on the pt integral. The fact that 
emissions below Qt are forbidden arises because 
the gluon not involved in producing the Higgs 
completely screens the colour charge of the fusing 
gluons if the wavelength of the emitted radiation 
is long enough, i.e. if pt < Qt- Now we see 
how this helps us solve our infra-red problem: as 
Qt — > so the screening gluon fails to screen and 
real emission off the fusing gluons cannot be sup- 
pressed. To see this argument through to its con- 
clusion we realise that multiple real emissions ex- 
ponentiate and so we can write the non-emission 
probability as 



exp 
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As Qt — > the exponent diverges and the 
non-emission probability vanishes faster than any 
power of Qt- In this way our integral over Qt be- 
comes 



-p^ L f(xi,Q T )f{x2, Qt) e" 



(9) 



which is finite. 

Now Eq. © is correct only so far as the lead- 
ing double logarithms. It is of considerable prac- 
tical importance to correctly include also the sin- 
gle logarithms. To do this we must re-instate the 
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running of a s and allow for the possibility that 
quarks can be emitted. Including this physics 
means we ought to use 



exp 



^Up\ a s (p 2 T ) 
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jf 1 A dz[zP 93 (z) + ^P 99 (z)]^ 



where A = 2pT/rriH, and P gg {z) and P qg {z) are 
the leading order DGLAP splitting functions. To 
correctly sum all single logarithms requires some 
care in that what we want is the distribution of 
gluons in Qt with no emission up to mu, and 
this is in fact [5] 



f(x,Q T ) 
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dlnQ T 



-S/2 
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The integral over Qt is therefore 



dQ 



'j-f(x 1 ,Q T )f(x 2 ,Q T ) 



Q 



which reduces to Eq. ([9]) in the double logarithmic 
approximation. 

Before we can go ahead and compute the cross- 
section we need to introduce the idea of gap sur- 
vival. The Sudakov factor has allowed us to en- 
sure that the exclusive nature of the final state is 
not spoilt by perturbative emission off the hard 
process. What about non-perturbative particle 
production? The protons can in principle inter- 
act quite apart from the perturbative process dis- 
cussed hitherto and this interaction could well 
lead to the production of additional particles. 
We need to account for the probability that such 
emission does not occur. Provided the hard pro- 
cess leading to the production of the Higgs occurs 
on a short enough timescale, we might suppose 
that the physics which generates extra particle 
production factorizes and that its effect can be 
accounted for via an overall factor multiplying the 
cross-section we have just calculated. This is the 
"gap survival factor". The gap survival, S 2 , is 
thus defined by 

da(p+H+p\no soft emission) = da(p+H+p)xS 2 
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Figure 4. CDF data on CEP dijet production. 
Figures from [7]. 



where da(p+H+p) is the differential cross-section 
computed above. The task is to estimate S 2 . 
Clearly this is not straightforward since we can- 
not utilize QCD perturbation theory. That said, 
data on a variety of processes observed at HERA, 
the Tevatron and the LHC can help us improve 
our understanding of "gap survival" and to date 
the HERA and Tevatron data do support the idea 
of gap survival. For the purposes of this talk we 
will presume to know the gap survival factor and 
that S 2 = 3% for CEP at the LHC, e.g. see [6J 
for an overview. 

Recently, the CDF collaboration has observed 
a 6cr excess of CEP of dijets at the Tevatron [7]. 
The agreement with the theory (as implemented 
in the ExHuME monte carlo generator [5]) is very 
good, as illustrated in Fig.QJ CDF also sees a sup- 
pression of quark jets in the exclusive region (high 
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Rjj), in accord with theoretical expectations. 

2. Higgs: SM and MSSM 

Fig. [5] shows how the cross-section for produc- 
ing a SM Higgs boson varies with Higgs mass (and 
for different gluon distribution functions). The 
cross-section is small and leads to low production 
rates. That said, a SM Higgs with mass above 
120 GeV should be observable in the WW* chan- 
nel with 300 fb _1 of data (which is around 3 years 
of high luminosity running) [119] . The gold-plated 
fully leptonic channel has very low backgrounds 
and has the advantage that one can still use the 
forward detectors to measure the mass. 

The bb channel is much more challenging. Trig- 
gering in this case would certainly benefit from 
having 220m detectors in place but even then one 
relies on optimistic scenarios for the production 
cross-section, detector acceptance and trigger effi- 
ciency. Nevertheless, it ought to be born in mind 
that CEP may be the only way to explore this 
channel at the LHC. 
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Figure 5. CEP cross-section for Standard Model 
Higgs producton. Figure from pQ. 



In contrast, the bb channel becomes much more 
exciting in certain MSSM scenarios. Rates are 
strongly enhanced at large tan/3 and small rriA, 



and the potential to measure the hbb Yukawa 
coupling is a strong selling point for CEP and 
a pre-requisite to determining any Higgs-boson 
coupling at the LHC (rather than just ratios of 
couplings) . Fig. [6] shows the region of parameter 
spac^l in which one could observe h — > bb using 
CEP [10] with different amounts of integrated lu- 
minosity. A similar pair of plots can be produced 
for H — > bb, see [TU]. Fig. [7] shows the result of 
an in-depth analysis of one particular point in the 
m^— tan/3 plane (tan j3 = 40 and vtia = 120 GeV) 
[T3] . The details of the two analyses can be found 
in |l|10|13j but the key point is that they are in 
general agreement 

The curves in Fig. [JJ correspond to different 
trigger scenarios. They also indicate the influence 
of pile-up and, in particular, the overlap back- 
ground (OLAP), in which the signal is faked by a 
coincidence of events, one that produces the cen- 
tral system (which fakes the Higgs decay) and one 
or more diffractive events that are able to produce 
protons in the forward detectors, e.g. a three- 
fold co-incidence of two single diffractive events 
with a pp — ► bbX event. Use of fast-timing detec- 
tors allows a significant reduction in the OLAP 
background, as the primary vertex can be pin- 
pointed to high accuracy. Improvements in the 
fast-timing could potentially eliminate the OLAP 
background completely and allow a 5c discov- 
ery with 3 years of high luminosity data taking 
(the mass peak is illustrated in the upper plot in 
Fig. EJ. 

3. Higgs: NMSSM 

To conclude, I would like to take a slightly 
more in-depth look at the possibilities for CEP of 
NMSSM Higgs bosons. More details can be found 
in [14]. The NMSSM is an extension of the MSSM 
that solves the /i-problem, and also the little hier- 
archy problem, by adding a gauge-singlet super- 
field S to the MSSM such that the \i term is now 
dynamical in origin, arising when the scalar mem- 
ber of S aquires a vev. The /i problem is solved 
since /x is no longer fundamental and therefore no 
longer naturally of order the GUT scale (as is the 
case if it is the only dimensionful parameter in 

3 In the M™ ax scenario with (j, = +200 GeV. 
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Figure 6. 5<r discovery contours (upper plot) and 
contours of 3er statistical significance (lower plot) 
for the h — > bb channel in CEP in the Ma - tan/3 
plane of the MSSM within the Mh,max benchmark 
scenario for different luminosity scenarios as de- 
scribed in the text [TD]. The values of the mass 
of the light CP-even Higgs boson, Mh, are indi- 
cated by contour lines. No pile-up background 
assumed. The dark shaded (blue) region corre- 
sponds to the parameter region that is excluded 
by the LEP Higgs boson searches [11112) . Figure 
from pp. 



Figure 7. Upper: Typical mass fit for the 
120 GeV MSSM h -> bb for 3 years of data tak- 
ing at 10 34 cm~ 2 s _1 after removing the over- 
lap background contribution completely with im- 
proved timing detectors. The significance is 5cr 
for these data. Lower: Significance of the mea- 
surement of the 120 GeV MSSM Higgs boson ver- 
sus luminosity, for two different combinations of 
muon (MU6, MU10) and fixed-jet-rate (FRT25, 
FRT10) triggers and with an improved (baseline) 
FP420 timing design (OLAP labels). Figure from 
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the superpotential) . The little heirarchy problem 
is also solved because a lighter Higgs is allowed, 
thereby taking the pressure off the stop mass. 
More specifically the lightest scalar Higgs can de- 
cay predominantly to two pseudo-scalar Higgses 
and the branching ratio to 6-quarks is correspond- 
ingly suppressed, thereby evading the 114 GeV 
bound from LEF0. Having a lighter Higgs means 
that the stop mass does not need to be so large, 
and that is preferred given the value of Mz- 

The Higgs sector of the NMSSM extends that 
of the MSSM by adding an extra pseudo-scalar 
Higgs and an extra scalar Higgs: crucially S is 
a gauge singlet and hence h —> aa can domi- 
nate with a light a (i.e. below the threshold for 
a — > bb). Freed of the heavy stop, it is most natu- 




draw attention to: 1. a light Higgs is preferred 
by the precision electroweak data (recall the best 
fit value is somewhat below 100 GeV); 2. a 100 
GeV Higgs with a reduced (10%) branching ratio 
to 6-quarks naturally accommodates the existing 
2.3cr LEP excess in e+e" -> Zbb |16ll7j . 



tan|8=3, M t 2 3 =100,200,300 GeV 




Figure 8. CEP of an NMSSM Higgs. 
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Figure 9. r c fr in units of 10~ 4 GeV 2 versus mh 
for tan/3 = 3. Point coding: (red) diamonds 
= 2m T < m a < 7.5GeV; (green) squares = 
7.5GeV < m a < 8.8GeV; (black) pluses = 
8.8GeV < m a < 2m b . Figure from [T4"] . 



ral to have a light Higgs with a reducing branch- 
ing ratio to &-quarks, as illustrated in Fig. [5J F 
and G are measures of fine-tuning, so the points 
in this scatter plot are supposed to represent most 
natural scenarios in the NMSSM. Our attention 
will focus on one such point, with m/, = 93 GeV 
and m a = 9.7 GeV with BR(h -> aa) = 92% 
and BR(a -> rr) = 81% 15 . The lightness of 
the pseudo-scalar a means that the h decays pre- 
dominantly to four taus. Should such a decay 
mode be dominant at the LHC, standard search 
strategies would fail and, as we shall see, CEP (as 
illustrated in Fig. |SJ) could provide the discovery 
channel. This "natural" scenario of the NMSSM 
has two additional bonus features that one might 

4 It drops to 86 GeV. 



To detect the four-tau decay of an NMSSM 
Higgs using CEP, we need first to trigger the 
event and to that end demand that at least one 
of the taus decays to produce a sufficiently high 
Pt muon. The muon then defines a vertex which 
can be used, in conjunction with the (picosecond) 
fast timing of the 420m detectors, to reject pile- 
up related backgrounds. The detailed analysis is 
outlined in [T3], here we shall just highlight the 
key features. Table 1 shows how the signal (CEP) 
and backgrounds (DPE, OLAP and QED) are af- 
fected by the cuts imposed. The top line of the 
table is the cross-section after imposing that there 
be at least one muon with px > 6 GeV, which is 
the nominal minimum value to trigger at level 
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Figure 10. The expected number of charged Figure 11. A typical a mass measurement, 
tracks reconstructed by the ATLAS inner detec- 
tor with pt > 6 GeV and \r]\ < 2.5. 



1 m ATLAS^ and the condition that both pro- 
tons be detected in the 420m detectors. There 
is also a loose cut on the invariant mass of the 
central system. Of the remaining cuts, I would 
like to single out the 'W ch = 4 or 6" cut. The 
charged track (iV c h) cut is noteworthy because it 
can be implemented at the highest LHC luminosi- 
ties: we cut on exactly 4 or 6 charged tracks that 
point back to the vertex defined by the muon. 
Pile-up events do add extra tracks (to both the 
signal and background), but they do not often 
coincide with the primary vertex (i.e. within a 
2.5 mm window) and do not spoil the effective- 
ness of this cut. The number of charged tracks 
in signal and background events is illustrated in 
Fig. QUI The ability to make such hard cuts on 
charged tracks could be of a much wider utility 
than this analysis (e.g. in defining a jet veto in 
Higgs plus two jets production). The four or six 
track event is then analysed in terms of its topol- 
ogy and the topology cut exploits the fact that the 
charged tracks originate from four taus in the sig- 
nal, which themselves originate from two heavily 
boosted pseudo-scalars. To avoid the effects of 
pile-up, the analysis is heavily track-based with 
almost no reliance on the calorimeter. 

5 It will turn out that a higher cut of 10 GeV is preferred 
in the subsequent analysis. 



Accurate measurement of the proton energies 
allows one to constrain the kinematics of the cen- 
tral system (in particular its invariant mass and 
mean rapidity are known). We can also extract 
the masses of the h and the a on an event-by- 
event basis. The mass of the h is straightforward 
of course (it is measured directly by the forward 
detectors) and a precision below 1 GeV can be 
obtained with just a handful of events. The mea- 
surement of the pseudo-scalar mass is more inter- 
esting and potentially very important. The pro- 
ton measurements fix p z and p x ^ y « for the cen- 
tral system. In addition, the tau pairs are highly 
boosted, which means they are collinear with 
their parent pseudo-scalars. That means that the 
four-momentum of each pseudo-scalar is approxi- 
mately proportional to the observed (track) four- 
momentum. The two unknown constants of pro- 
portionality (i.e. the missing energy fractions) are 
overconstrained, since we have three equations 
from the proton detectors. The result is that we 
can solve for the pseudo-scalar masses, with four 
measurements per event. Fig. [TT] shows a typical 
distribution of a masses based on 180 fb _1 of data 
collected at 3 x 10 33 cm^V 1 . 

In Table 2 I show the bottom line numbers for 
three different trigger scenarios and three differ- 
ent instantaneous luminosities. The key point is 
that, although we can expect only a handful of 
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Table 1 

The table of cross-sections for the signal and backgrounds. All cross-sections are in femtobarns. The 
overlap (OLAP) background is computed at a luminosity of 10 34 cm~ 2 s _1 . 





CEP 


DPE 


OLAP 


QED 


Cut 


H 


bb 


99 


bb 


bb 


4r 


2r 21 


Pro, fi, 6) M 


0.442 


25.14 


1.51xl0 3 


1.29xl0 3 


1.74xlO B 


0.014 


0.467 


N ch = 4 or 6 


0.226 


1.59 


28.84 


1.58xl0 2 


1.44xl0 4 


0.003 


0.056 


Q h = 


0.198 


0.207 


3.77 


18.69 


1.29xl0 3 


5xl0~ 4 


0.010 


Topology 


0.143 


0.036 


0.432 


0.209 


1.84 




<0.001 


Py, isolation 


0.083 


0.001 


0.008 


0.003 


0.082 






Pr 


0.071 


5xl0~ 4 


0.004 


0.002 


0.007 






m a > 2m T 


0.066 


2xl0~ 4 


0.001 


0.001 


0.005 







signal events, the background is under control. 
Remember, we need only a few events in order 
to extract the masses of both the a and the h. 
The statistical significance of any discovery is es- 
timated in Fig. 1121 The lower plot might pertain 
if the signal rate were doubled (recall the theoret- 
ical uncertainty permits it) or if data from AT- 
LAS and CMS were combined or if other leptonic 
triggers are used pQ . 

4. Concluding remarks 

Central exclusive production is a very attrac- 
tive prospect for the LHC. A very broad pro- 
gramme of physics can be pursued for very little 
additional expenditure. Measurements from CDF 
at the Tevatron are very encouraging and support 
the validity of the theoretical calculations: the 
theory is probably not too far off the mark. More- 
oever, one of the highlights of the past couple of 
years has been the demonstration that high lumi- 
nosity backgrounds can be brought under control. 
This talk has focussed only upon Higgs physics 
and has placed particular emphasis on the possi- 
bility that CEP could be the only way one could 
observe at the LHC a Higgs boson that decays 
predominantly to four taus: something that could 
be fairly generic feature of SUSY models with an 
enlarged Higgs sector (such as the NMSSM). 
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Table 2 

Expected number of signal (S) and background 
(B) events for the three trigger scenarios assum- 
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taneous luminosity of lxlO 33 cm~ 2 s _1 (row 1), 
5xl0 33 cm~ 2 s _1 (row 2) and 10 xlO 33 cm~ 2 s" 1 
(row 3). 
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